Solitary coherent structures in viscoelastic shear flow: computation and 

mechanism 



K. Arun Kumar and Michael D. Graham 
Department of Chemical Engineering and Rheology Research Center, University of Wisconsin-Madison, Madison, WI 

53706-1691. 
(February 2, 2008) 

ABSTRACT 

Starting from stationary bifurcations in Couette-Dean flow, we compute nontrivial stationary solutions in in- 
ertialess viscoelastic circular Couette flow. These solutions are strongly localized vortex pairs, exist at arbitrarily 
large wavelengths, and show hysteresis in the Weissenberg number, similar to experimentally observed "diwhirl" 
patterns. Based on the computed velocity and stress fields, we elucidate a heuristic, fully nonlinear mechanism 
for these flows. We propose that these localized, fully nonlinear structures comprise fundamental building blocks 
for complex spatiotemporal dynamics in the flow of elastic liquids. 

PACS. 47.54.+r Pattern formation - 47.50.+d Non-Newtonian fluid flows - 83. 50. Ax Steady shear flows - 83.85.Pt 
Flow computation. 



Flow instabilities and nonlinear dynamics have long been recognized to occur in flows of viscoelastic poly- 
mer melts and solutions jj]-^. An important breakthrough, which has led to increasing recent attention to 
these phenomena, was made by Larson, Shaqfeh and Muller 0, who discovered that circular Couette flow of a 
viscoelastic liquid undergoes an instability loosely analogous to the classical Taylor- Couette instability of New- 
tonian liquids, but driven solely by elasticity - the instability is present at zero Taylor (Reynolds) number. They 
also showed that a linear stability analysis of a simple fluid model predicts instability and elucidated the basic 
elasticity-driven mechanism of the instability. In particular, although many researchers had studied the effects 
of viscoelasticity on the Newtonian (inertial) Taylor-Couette instability, both experimentally [||-^), and theo- 
retically |lC| , pd] ,|5|-[7t^2| , these workers were the first to demonstrate an inertialess, purely elastic mechanism for 
instability in a viscometric flow. More recent observations have revealed a wealth of interesting dynamics in this 
flow as well as other simple flows ||[l3]-[ll| ; one set of experimental observations of particular interest was made 
by Groisman and Steinberg [ |l3[ ], who found long- wavelength, stationary axisymmetric vortex pair structures in 
inertialess viscoelastic flow in the circular Couette geometry. There are two interesting aspects to these obser- 
vations: (1) isothermal linear stability analysis in this geometry never predicts bifurcation of stationary states 
(in contrast to the classical Taylor-Couette case) and (2) the observations suggest that these vortex pairs, which 
Groisman and Steinberg dubbed "diwhirls" , can exist in isolation - there does not seem to be a selected axial 
wavelength for this pattern. These considerations motivate the present computational study, which addresses 
the following questions: (1) Do isolated branches of stationary solutions exist in a simple model of a viscoelastic 
fluid in the circular Couette geometry? (2) If so, what are the spatial structures of these? Are they localized? 
(3) A nonlinear self-sustaining mechanism must be present for such patterns to exist. Can the computations 
help us elucidate it? 

We address these questions by fully nonlinear computations of the branching behavior of an inertialess isother- 
mal FENE dumbbell fluid in the circular Couette geometry. Introduction of an azimuthal body force as an 
additional parameter allows the possibility of accessing steady solution branches that would be isolated in the 
parameter space of Couette flow. We find in short that the answers to all of the above questions is yes: isolated 
solution branches consisting of spatially localized coherent structures have been computed for the first time in a 
purely clastic flow, in very good agreement with the experimental observations, and a heuristic, fully nonlinear 
mechanism has been elucidated. 



We consider the flow of an inertialess polymer solution in the annulus between two concentric cylinders. The 
inner cylinder has radius R\ and the outer cylinder has radius i?2- The flow is assumed to periodic in the axial 
direction, and we denote the period, nondimensionalized with the gap width R 2 — R\, by L. The fluid has a 
relaxation time A; the polymer and solvent contributions to the viscosity are denoted respectively by rj p and T} a , 
with the ratio T) s /r) p denoted by S. The flow is created by a combination of the motion of the inner cylinder at 
a velocity £lRi and by the application of an azimuthal pressure gradient Kg = dP/d6. The equations governing 
the flow are the momentum and mass conservation equations, and the FENE-P constitutive equation |lq ], which 
models polymers as beads connected by finitely extensible springs. These may be written in dimensionless form 
as 

V ■ r - Vp + We 9 SV 2 v = 0, (1) 
V • v = 0, (2) 

We e (°£- {a . V.}t - {« • V.}) + ° Q)/6) - l) = 0, (3) 

where v is the velocity, p is the pressure, a is the ensemble average of the conformation tensor, and r = 
a/(l — tr(a)/b) — I is the polymer stress tensor, with \fb being a dimensionless measure of the maximum 
extensibility of the dumbbells. The Weissenberg number, Wee is the product of the polymer relaxation time 
and a characteristic shear rate, which we take to be the shear rate at the outer cylinder for an Oldroyd-B 
fluid (1/6 = 0) flowing through the geometry. Other important parameters are the dimensionless gap width, 
e = (i? 2 — i?i)/i?2, and 5, which measures the relative importance of the pressure gradient as the driving force 
for the flow, given by 

s = -K g e 2 R 2 /(2 (r, p + r, 8 )) 

(l-e)R 2 n-K e e 2 R 2 /(2(r 1p + r ]s )y {) 

so that S = is circular Couette flow and 5 = 1 is Dean flow (pressure driven flow in a curved channel). The 
velocity satisfies no slip boundary conditions on the walls of the cylinder. We compute steady, axisymmct- 
ric solutions to the governing equations using a spectral element scheme JU7 with Galerkin weighting on the 
conservation equations and streamline-upwind/Petrov-Galerkin weighting J^ji| on the constitutive equations. 
We use a branch switching technique to compute starting points on non-trivial branches, and pseudo-arclength 
continuation to trace out solution curves in parameter space [ p0[ . The linear systems in the Newton iteration 
are solved using GMRES |21| , with a modified dual threshold preconditioner f2^,^j| . 

We restrict our attention to parameters which are close to those used in the experiments by Groisman and 
Steinberg ]13|| , Specifically, we choose parameter values S — 1.2, and e = 0.2. The Reynolds numbers in their 

experiments were ~ 1, so our neglect of inertia corresponds well to their experiments. In our computations, we 
find that the most interesting behavior occurs only for large values of 6, so for the results we present below, we 
choose a value of 1830. 

In the Couette flow case (5 = 0), there is no bifurcation from the viscometric base state to a stationary 
nontrivial state - if they exist, any branches of nontrivial stationary states in this case are isolated. So to 
pursue such solutions, an indirect approach must be taken. We use the fact that in pure pressure-driven flow 
(5 = 1), stationary bifurcations do exist - by continuity these bifurcations also exist for 6 < 1 and we search for 
a route to the Couette case by exploring the evolution of these bifurcating branches. Specifically, we traced out 
a stationary branch originating at S — 0.576, L = 2.71, Wee — 25.15, then increased L to 3.08 while holding the 
values of the other parameters fixed, then reduced 5 downward, finding that this solution branch persisted at 
5 = 0. We have thereby shown that there is indeed an isolated branch of nontrivial stationary patterns in the 
circular Couette flow geometry. If we now decrease L, we find a turning point at L = 2.94, so these solutions do 
not exist for wavelengths shorter than this. On the other hand, when we increase L, we find that the solutions 
exist at arbitrarily large wavelength. Figure 1 shows the contours of the streamfunction and aee at L = 110.89 
( i.e., the axial wavelength is more than 100 times times the gap width!) The core of the pattern (which has 
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been magnified for clarity) is a region of very strong inflow, surrounded by regions of much weaker outflow. The 
aee field shows an even more intense localization at the core and requires the use of a very fine spectral element 
mesh to resolve. The strong localization at the core and the asymmetry between inflow and outflow are features 
shared by the experimentally observed diwhirls fL3| ], which indicates that we are capturing the physics behind 
the structures they observe. Away from the core, the flow field is almost pure circular Couette flow. In obtaining 
solutions at very large values of L, we exploited the lack of axial variation away from the center and the fact 
that the characteristic length of the localization near the core does not change at large wavelengths by using 
the same mesh as at smaller values of L close to the center, and simply increasing the lengths of the spectral 
elements bordering the axial edges. 

Figure 2 shows the results of continuation in Weg for three different values of L. In each case, we see a turning 
point in Weg. It is located at Weg = 24.97 for L = 3.08, at We e = 23.37 for L = 4.72, and at Weg = 23.55 
for L = 9.11. The point at which the diwhirls lose existence is much lower than the critical Weg for linear 
instability at that wavelength, so that the overall bifurcation structure shows a hysteretic character, consistent 
with experiment. 

A natural question regarding these new flow states is their stability. We have performed linear stability 
analyses with respect to axisymmetric perturbations, for selected parameter values. These results will be detailed 
elsewhere p3[ , but the main result is simple and interesting: there is an unstable complex conjugate pair of 
eigenvalues, corresponding to instability with respect to oscillatory disturbances. The destabilizing disturbance, 
however, has significant amplitude only near the ends of the domain, where the flow is essentially Couette flow 
(which is linearly unstable at the parameter values chosen) and vanishes in the core region of the diwhirl. So the 
picture that emerges is this: the diwhirl pattern is robust, and it coexists with, but is spatially distinct from, 
the oscillatory finite wavelength pattern arising from the linear instability. 

By estimating the work done by a fluid element in a single flow cycle (as seen in a Lagrangian reference 
frame), Groisman and Steinberg E4| argued that long wavelength stationary structures which exhibited significant 
asymmetry between inflow and outflow were possible in Couette flow. While this argument shows that diwhirl 
structures are physically plausible, it does not explain the mechanism by which they are sustained. Having the 
complete velocity and stress field available to us, we propose a more complete mechanism. Figure 3 shows a vector 
plot of v near the outer cylinder at the center of the diwhirl structure. We see that the vg field at the diwhirl 
center is locally parabolic near the outer cylinder, similar to the velocity field in the upper half of the channel in 
Dean flow, which displays a stationary vortex instability due to the corresponding unstable stratification of the 
azimuthal normal stress |2q] . We therefore propose the following mechanism for the instability: a local, finite 
amplitude perturbation near the outer cylinder creates a locally parabolic velocity profile. This velocity profile 
results in an unstable stratification of hoop stress, just as in Dean flow, which drives inward radial motion. The 
fluid accelerates azimuthally as it moves radially inward, due to the base state velocity gradient, thus maintaining 
the increase in vg, resulting in a self-sustaining process (figure 4). 

To summarize, we have computed stationary nontrivial solutions in circular Couette flow using a FENE-P 
model by numerical continuation from stationary bifurcations in Couette-Dean flow. These solutions show very 
strong localization, exist only for large values of b and large wavelengths, and show a hysteretic character in Weg. 
The self-sustaining mechanism is related to the mechanism of instability in viscoelastic Dean flow, and arises 
from a finite amplitude perturbation giving rise to a locally parabolic profile of the azimuthal velocity near the 
upper wall. The computed flow structures are very similar to the experimentally observed diwhirl patterns p"3] | . 
Along with the solutions arising from the linear instability of the circular Couette flow base state, we propose 
that these solutions form building blocks for spatiotemporal dynamics in the flow of elastic liquids. 
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Khomami, Prof. R. Sureshkumar, and their research associates for helpful discussions. 



[1] C. J. S. Petrie and M. M. Denn. AIChE J., 22(2):209-236, 1976. 
[2] R. G. Larson. Rheol. Acta, 31:213-263, 1992. 



3 



E. S. G. Shaqfeh. Ann. Rev. Fluid Mech., 28:129-185, 1996. 

R. G. Larson, E .S. G. Shaqfeh, and S. J. Muller. J. Fluid Mech., 218:573-600, 1990. 

H. Rubin and C. Elata. Phys. Fluids, 9:1929-1933, 1966. 

M. M. Denn and J. J. Roisman. AIChE J, 15:454-459, 1969. 

Z.-S. Sun and M. M. Denn. AIChE J, 18:1010-1015, 1972. 

H. Giesekus. Prog. Heat Mass Transfer, 5:187-193, 1972. 

G. S. Beavers and D. D. Joseph. Phys. Fluids, 17:650-651, 1974. 

S. K. Datta. Phys. Fluids, 7:1915-1919, 1964. 

R. H. Thomas and K. Walters. J. Fluid Mech., 18:650-651, 1964. 

M. Lange and B. Eckhardt. Preprint, 2000. Although these authors denote the solutions that they find as "diwhirls", 
their investigation is confined to regimes where Re 3> 1 and Wee /Re <C 1, where Re is the Reynolds number. In 
these parameter regimes, fluid inertia plays a significant role. Thus, the patterns they simulate cannot be classified 
as being driven by elasticity alone. 

A. Groisman and V. Steinberg Phys. Rev. Lett, 78(8): 1460-1463, 1997. 

B. M. Baumert and S. J. Muller. J. Non-Newtonian Fluid Mech., 83(l-2):33-69, 1999. 
A. Groisman and V. Steinberg. Nature, 405:53-55, 2000. 

R. B. Bird, C. F. Curtiss, R. C. Armstrong, and O. Hassager. Dynamics of polymeric liquids, volume 2. Wiley, New 
York, 2nd edition, 1987. 

Y. Maday and A. T. Patera. In State of the art surveys on computational mechanics, pages 71-143. ASME, 1989. 

A. N. Brooks and T. J. R. Hughes. Comp. Methods Appl. Mech. Eng., 32:199-259, 1982. 

J. M. Marchal and M. J. Crochet. J. Non-Newtonian Fluid Mech., 26:77-114, 1987. 

R. Seydel. Practical bifurcation and stability analysis. Springer- Verlag, New York, 1994. 

Y. Saad and M. H. Schultz. SI AM J. Sci. Stat Comput., 7(3):856-869, 1986. 

Y. Saad. Iterative methods for sparse linear systems. PWS Publishing Company, Boston, 1996. 

K. A. Kumar and M. D. Graham. In preparation. 

A. Groisman and V. Steinberg. Phys. Fluids, 10(10):2451-2463, 1998. 

Y. L. Joo and E. S. G. Shaqfeh. Phys. Fluids A, 3(7):1691-1694, 1991. 



4 



FIGURE CAPTIONS FOR KUMAR AND GRAHAM 



Figure 1: Density plot of agg (white is large tension, black small) and contour plot of the streamfunction at 
L = 110.89 {Weg = 23.52, b = 1830, S = 1.2, and e = 0.2). For clarity, most of the flow domain is not shown. 
Note the very strong localization of age near the center. Away from the core, the structure is pure circular 
Couette flow. 

Figure 2: Diwhirl solution amplitudes as functions of Wee and L. Note that the curves at L = 9.11 and L = 4.72 
are very close together, while both curves are well separated from the curve at L = 3.08 (b = 1830, S — 1.2, and 
e = 0.2) 

Figure 3:Vector plot of v near the outer cylinder at the center of the diwhirl structure (oblique arrows) and 
the base state (straight arrows) . The length of the arrows is proportional to the magnitude of the velocity. The 
axial velocity is identically zero in the base state, and is zero by symmetry at the center of the diwhirl. 

Figure 4: Nonlinear self-sustaining mechanism for the diwhirl patterns. 
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